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Abstract 

^"^^ ' Applying a unifying Lax operator approach to statistical systems a new class of integrable vertex 

' models based on quantum algebras is proposed, which exhibits a rich variety for generic q, q roots 

■ of unity and q 1. Exact solutions arc formulated through algebraic Bethe ansatz and a novel 



possibility of hybrid vertex models is introduced. 



^ • PACS numbers 02.30.Ik, 02.20.Uw, 05.20.-y 05.50+q, 03.65.Fd 
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CJ ■ Introduction: 

D-dimensional quantum systems are known to be related to (l+D)-dimensional classical statistical 
^ ' models, which is true naturally also in D = 1, where one finds an exclusive class of models, known as 
■ integrable systems, allowing exact solutions. Celebrated example of such relation is that between the 
XYZ quantum spin-^ chain and the 8- vertex statistical model and similarly between the XX Z chain 
and the 6- vertex model [Q]. Hamiltonian Hs of the integrable quantum spin chain is given through 
its transfer matrix as lnT{u) = I + uHg + O(u^), while the partition function Z of the related vertex 
model is constructed from t{u) as Z = tr{T{u)^'^). Moreover, both these models share the same 
quantum i?-matrix and have the same representation for the transfer matrix, commutativity of which: 
[t[u),t{v)] = guarantees their integrability. 

It is therefore rather surprising to note that, inspite of such deep connection between these two in- 
tegrable systems, their starting formulation conventionally follows two different routes. Quantum sys- 
tems usually are defined by their Lax operators Lai (u) , which satisfy the quantum Yang-Baxter equa- 
tion (YBE) Rab{u — v) Lai{u) Lbi{v) = Lbi{v)Lai{u)Rab{u — v), together with its associated i?-matrix. A 
vertex model on the other hand is described by its Boltzmann weights given generally through the ele- 
ments of the i?-matrix alone, which solves the YBE Rab{u — v)Rai{u)Rbi{v) = Rbi{v)Rai{u)Rab{u — v). 
However such a difference in their approaches, reason of which seems to be rather historical, puts 
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certain restrictions on the 2-dimensional vertex models by assuming their vertical (v) and the hori- 
zontal (h) links, which are related to the auxiliary and quantum spaces respectively, to be equivalent. 
As a consequence while a rich variety of integrable quantum systems with wide range of interactions 
involving spin, fermionic, bosonic and canonical variables does exist, the integrable vertex models 
are confined mostly to those quantum models that exhibit regularity property expressed through the 
permutation operator: Lai{0) = Pal and hence correspond to local Hamiltonians with nearest neighbor 
(NN) interactions. Therefore the well known examples of the integrable vertex models, apart from 
those mentioned above, appear to be limited mainly to the models like the 5-vertex model Q, 6- vertex 
model in external fields Q, 19- vertex model connected with the Babujian-Takhtajan spin-1 chain |Q] 
and the vertex models equivalent to the Hubbard model ||5|, supersymmetric t-J model [^, Bariev 
chain Q etc. all exhibiting only NN interactions. 

The basic idea of the present letter however is to exploit fully the equivalence between statistical and 
quantum systems and construct new class of integrable vertex models by applying a unifying scheme 
designed originally for quantum models [§. In the original scheme an ancestor model was proposed 
for generating integrable quantum systems as its various descendant realisations. For describing our 
vertex models we start in analogy also with the generalised Lax operator 

^ / c+g^'+" + qg-(^'+") 2sina5- 

V 2sinQ5+ + C2 9^'"'' 

linked with the underlying quantum algebra 

[S^S^]=±S^, [S~^,S-] = {M+[2S% + M-[[2S%), [M±,-]=0. (2) 

Here fxL = [[xll„ = '^"'^^"^^ and the central elements are related to the other set of such 

elements appearing in the L-operator as = ^^^/^l{c~^cZ ± c^c^). It is important to notice that 
(§) is a q-deformed quadratic algebra, which generalises both q-spin and q-boson algebras and in fact 
follows from the quantum YBE representing integrability condition. We define the Boltzmann weights 
(BW) of our vertex models not by the i?-matrix as conventional, but through the elements of the 
Lax operator: ^^^{u) = u}a,j-b,k{u) by using matrix representations of the general algebra (§) in (||). 
These generalised BW generate a unified vertex model, which through possible reductions yields new 
series of vertex models linked with different underlying algebras, their representations and choices of 
the central elements. Prominent examples of such integrable statistical systems are a rich collection 
of vertex models linked to q-spin and q-boson with generic q, q roots of unity and q ^ 1. In all these 
models the h and v links , contrary to the usual approach, may become inequivalent and independent 
at every vertex point and since we consider here 2x2 Lax operators, the h links admit only 2 values: 
a,b £ [+)—]• The v links on the other hand have richer possibilities with j,k £ [^,D], depending 
on dimension D of the matrix-representation of the q-algebras (see Fig. 1). The familiar ice-rule is 
generalised here as the 'colour' conservation a + j = b + k for determining nonzero BW. The crucial 
partition function of the models however is given as usual hy Z = J2configHa,b,j,k^a,j;b,kiu). 

An important point to note is that unlike traditional approach the Lax operators related to such 
vertex models do not coincide with their i?-matrix, do not comply with the regularity condition and do 
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not correspond in general to quantum Hamiltonians with NN interactions. Moreover since our vertex 
models belonging to the same class have the same i?-matrix, we can generate another rich series of 
integrable models, namely hybrid vertex models by combining any number of them in a row (see Fig. 
!)• 

The eigenvalue solution of the transfer matrix needed for constructing the partition function for 
all these vertex models can also be found exactly through the algebraic Bethe ansatz in a unifying 
way. 

Unified vertex model: 

In accordance with our primary goal we discover an explicit matrix representation for the basic 
operators S^, S^: 

< s,m\S^\m,s >= m6m,m, < s,m\S'^\m, s >= f^{'m)6m±i,m, (3) 

with /+(m) = f-{m + I) = {k + [s - m]g{M+[s + m + l]g + M-[[s + m + l]]q))^ having additional 
parameters k, s. It may be checked that (^) indeed gives an exact representation of the general q- 
deformed algebra (^) for arbitrary values of the central elements M^. Therefore the BW may be 
constructed from the matrix representation of the generalised Lax operator (|^) by using in the 
form 

^±,,;±j(n) = 4e^°("±"^) + cie-^"^'^^™), a;+,,-;_,,-_i = = 2/+(m) sin a, (4) 

with m = s + 1 — j, j = 1,2, . . . , D. The BW parameterised as would generate now a unified 
(4D — 2)-vertex model, representing a new series with arbitrary parameters c^, s and k. These models 
and naturally all others obtained below through various reductions are integrable statistical models 
and share the same i?-matrix, which is given by that of the well known 6- vertex model . 

Note that though in general the dimension D of the matrices ^ is infinite, it may get truncated 
through possible appearance of zero-normed states. To analyse this important effect we observe that 
since [0]q = 0, one gets f^im = s) = for k = 0, recovering the familiar 'vacuum' state: S~^\s, s >= 0. 
However due to the presence of the second term one gets here f^{m = —s) = {[2s + l]q(M"''[0]q + 
M~[[0]]g))2 ^ and unlike the spin representation we have S~\m,s >^ for all m. This creates 
therefore an infinite tower of states by the action of the lowering operator S~ , as typical for the 
bosonic representation. This also signals the fact that algebra (|2|) includes both q-spin and q-boson 
and therefore their representations must be derivable from (^) as particular cases. 
q-spin vertex model: 

It is straightforward to check that for = 1,M~ = 0, our unifying algebra (|2|) reduces to the 
well known Uq{su{2)) quantum spin algebra ||9| and at the same time @ reproduces the known q-spin 
representation. Therefore the corresponding BW may be obtained from for a consistent choice 

= = =F^i as 

'^±,j;±,j{u) = [u± m]q, uj+j.-j-i = = /+(^^P™)(m), m = s + 1 - j 

with ft^'^^^^^\m) = ([s =F m]q[s it m + l]g)2 . In this case the truncation S^\m = is, s >= typical 
for spin models and hence the familiar D = 2s + 1 dimensional representation naturally arise, which 
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produces therefore a series of q-spin (8s + 2)-vertex models. The 6-vertex model is clearly recovered 
at s = ^, while s = 1, |, 2, |, . . . yield new 10, 14, 18, 22, . . .-vertex models (Fig. 1 a,b)). 

The quantum systems related to such statistical models may be represented in general by inter- 
acting q-spins with nonlocal interactions. In particular, since the well known sine-Gordon model is a 
realisation of the q-spin |10|, the vertex models constructed with nonzero k and having infinite D will 
be related to the quantum integrable lattice sine-Gordon model |11]. 

q-boson vertex model: 

We find that q-bosonic algebra jl^ can also be derived as a subalgebra of (^) for the complementary 
choice M"*" = sina,M~ = cos a by denoting = pA, S~ = pA^ , = —N, p = (cot a) 2. The 
same choice derives therefore the matrix representation for the q-boson directly from (^) with the 
assumption that k = s = and n = —m, yielding 

fo^''"'\n) = ([1 + nUi-n - 1]],)^ = -^[1 + n] J, ff'^'-\n) = f^^''-\n - 1) = -^N J. 

Consequently for a consistent solution cjl = 1, = =pze^*° we can derive from (^) the BW as 

UJ±J;±A^) = ie±*°'^[^zT(i + 0-l)]„ ^=\(^ + ^)^ 

= ^-,-i;+, = /o'^^''°^^(j-l) = ^b--4. (5) 

It is obvious that apart from the vacuum state |0 > with fQ^'^^°^\o) = "^[0]^2 = we can 
have no other zero-normed states and the q-bosonic representation like the standard boson is semi- 
infinite with D = n + 1. The integrable (4n -|- 2)-vertex model linked to the q-boson (Fig. Ic)) that 
we construct using (^) would therefore be related to the lattice version of the quantum derivative 
nonlinear Schrodinger model (DNLS), which exhibits a q-bosonic realisation []T3|. 



Vertex models with q roots of unity: 

An excellent possibility for regulating the dimension of the matrix representation opens up when 
q = is chosen as solutions of = ±1 with parameter a taking discrete values Oa = t^^jCi = 

1,2, ... ,p — 1 [14|. Note however that when some values of a becomes a factor of p one faces a 

p 

situation with qa = ±1. Therefore to avoid such complicacies we suppose p to be prime in our present 
discussion. For further analysis we focus on the action of S" assuming k = in (^): S~\m = —s, s >= 
([s -|- s -|- l]q(M"'"[s — s]q + M^[[s — and observe that due to [p]^ = sin OaP = 0, unlike generic 

q we can get now S~ \ — s,s >= at s = p — {s + 1), which reduces matrix (^) to a finite dimensional 
representation. Therefore the BW obtained from (^) for this case would generate another series of 
unified K- vertex model having finite K = 4p — 2 configurations at every vertex point. Moreover, since 
for a fixed p there can he p — 1 different Oa, each of these discrete values describes a different set of 
BW and hence a novel model. 

Consequently at particular reductions as analysed above, we obtain the corresponding series of new 
vertex models linked with q-spin or q-boson, but now having finite configuration space determined by 
p. Thus for the q-spin with fixed p, < p < 2s -|- 1, in place of a (8s -|- 2)-vertex model for generic q, 
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one obtains p — 1 number of different (4p — 2)-vertex models and the related representation including 
the case p > 2s + 1 become more involved |jl^ . The corresponding BW defining these models should 
however be given by their same generic form, though using discrete values. As for example in case 
of s = I with = —1, instead of a 22- vertex model one obtains 4 different 18- vertex models for 
distinct values of = vr#,a = 1,2,3,4. Noticeably, as a quantum model the q-spin with q roots of 



unity are realised as the restricted sine-Gordon model |15|. 

The situation becomes more interesting when applied to the q-boson with finite p, since now 
together with the standard vacuum we get also A^\n = | — 1 >= ■^[|]g2|f >= 0, yielding finite 
(| X |) matrix representations for the q-bosonic operators A,A^. As a result we obtain an intriguing 
series of {2p — 2)-vertex models with BW described by the same form (^) as for the generic q-bosonic 
case, but with different possible parameter values q = e*"", o = 1, 2, ... ,p — 1. The quantum model 
corresponding to such q-boson vertex models can be realised as the restricted DNLS model, which 
supports finite quasi-particle bound states 
Rational class of vertex models : 

At g — > l(a — > 0) on the other hand, the associated i?-matrix goes to its known rational limit 
and the underlying algebra becomes undeformed one with m^, reducing at the same time 

the unified model to its rational form. Consequently, taking carefully the limits we may construct 
in a similar way the corresponding set of vertex models belonging to the rational class. Not going 
into details we mention only that the BW of these vertex models can be obtained from the limiting 
values of (Q) yielding f^{m) — > ((s — m){rn^ {m -|- s -|- 1) -|- m~))2 . It is easy to check that the BW 
for the vertex models related to the undeformed spin as well as the standard boson correspond to 
the particular values of the central elements: m"^ = l,?7i~ = and m"*" = 0, m~ = 1, respectively. 
Remarkably, the spin vertex model constructed in this way coincides with the similar higher s model 
obtained earlier through fusion method whereas the bosonic- vertex model apparently is a new 
model, linked to a quantum integrable lattice NLS model p7[| . 
Hybrid vertex models: 

In constructing our vertex models we have flatly assumed that in any model the same BW must 
be defined at every vertex point. An immediate generalisation is therefore possible by relaxing this 
condition and considering the central elements as well as the spin parameters s appearing in (^) 
to be different at different sites. As we have already stressed, vertex models obtained as various 
reductions of the same integrable unified model belong to the same class sharing the same i?-matrix. 
Thus the q-spin and q-boson vertex models are members of the trigonometric class, while the normal 
spin and boson models belong to the rational class. Based on this fact therefore we can construct a 
rich collection of hybrid models by combining different vertex models of the same class and inserting 
their defining BW along the vertex points Z = 1, 2, . . . , in a row, in any but fixed manner. Due to 
the association with the same i?-matrix the integrability of such statistical models would be naturally 
preserved. 

Thus for example an alternate insertion of 10 and 6 vertex models results to a hybrid model, which 
is related to the known quantum model [|l^] involving spin-1 and spin-^ operators with next-NN 
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interactions. More exotic hybrid models can be formed by arranging the BW for the q-spin and q- 
boson vertex models, alternatively or in any other way at different vertex points (Fig. 1). Similarly one 
can construct a spin-boson hybrid vertex model by combining their individual vertex models, which 
would correspond to a quantum chain of interacting spins and bosons involving next-NN couplings. 

Unified solution: 

The construction of the unified vertex model through the generalised Lax operator suggests also a 
scheme for exactly solving the eigenvalue problem for the transfer matrix. Since the partition functions 
in turn can be determined from the knowledge of these eigenvalues, all vertex models obtained as 
particular cases and linked to (un-)deformed spin or (un-)deformed boson can also be solved in a unified 
way. There is a well formulated algebraic Bethe ansatz method for exactly solving the eigenvalue 
problem of the transfer matrix: t{u) = trh{Ylj^ Li{u)), when the associated Lax operator as well as 
the i?-matrix are given |19[ |. Therefore, since we have defined the BW through matrix representations 
of the Lax operator and the i?-matrix in our case is given by that of the well known 6-vertex model, 
we can derive the exact eigenvalues for the transfer matrix of our models as 



n n r _|_ 1 1 

A(n) =<i^+,i(n)n5(^fc-n)+u;V,i(^.)n5(u-nfc), g{u) = i^^^. (6) 

k k 



with all possible solutions of {uk} to be determined from the Bethe equations 

f^;^^i^\\f^ lu,-u, + l] (7) 

ya;_,i;_,i(nz) J [ui - Uk- l\q 

By analysing the structure of these equations we conclude that, the factors involving BW in both of 
them come from the action of the Lax operator on the pseudovacuum, which is chosen as the direct 
product of the highest weight states with j = 1 i.e. \m = s >. The rest of the factors on the other 
hand are originated from the i?-matrix elements, which arise during diagonalisation of the transfer 
matrix due to the use of the quantum YBA. Therefore it is crucial to note that, the only part given by 
BW is model-dependent and defined for the vertex models by the diagonal entries in (^) with j = 1, 
while the remaining parts contributed by the i?-matrix are the same for all our models from the same 
class. Consequently the exact solutions for all models constructed here can be found in a systematic 
way from (P) and (0) by using corresponding reductions of the unified model (^). 

The total number of solutions {A^(n)} for the eigenvalues (^) should be , which coincides with 
the number of possible eigenstates and gives the dimension of the vector space on which the transfer 
matrix acts. The partition function of the vertex models may therefore be given at the thermodynamic 
limit by Z = limM,N^ootrv{T^\u)) = limM,N^ooJ2^=i i^)- this important limit, the Bethe 
equations (0) turn into an integral equation V{u) = 2t: p{u) + J p{v)K(v,u)dv, with known kernel 
of the 6-vertex model |20|. Interestingly, all information about a particular model is encoded in the 



driving term only, which is expressed through = re*'^^"-' as V{u) = P' (u) and therefore 

knowing the explicit form of BM one can derive easily the equations for individual models. 
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For extracting the solutions of the hybrid vertex models however the BW dependent parts in the 
above equations should be slightly modified by generalising the factors inhomogeneously as Ylf^{{uj^^\.^ il^))^*^ i 
where Nf^ is the number of vertices of type /? appearing in a row with the constraint N = J2f3 ^p- 

Detail investigation of individual models and identification of their most probable states are im- 
portant problems to be pursued. 
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Figure 1: Integrable vertex models with horizontal (h) links taking 2 values, while the vertical (v) ones 
may have D possible values, a) 6-vertex b) q-spin vertex and c) q-boson vertex models. Combining 
a,b,c) an integrable hybrid model may be formed, = 1 gives D = p inh) and c) 
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